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1. Introduction 



When it was first introduced in its modern form flj, the question "Does unitarity and scale 
invariance imply conformal invariance?" was mostly of academic interest. Recent work [2jp) 
showed that scale-invariant theories display renormalization group (RG) flow recurrent be- 
haviors and have novel implications for beyond the standard model phenomenology [4]|^] 
Thus, the existence of scale-invariant theories has deep consequences, especially with respect 
to the intuitive understanding of RG flows as the integrating out of degrees of freedom, 
and the c-theorem. "Does unitarity and scale invariance imply conformal invariance?" is 
therefore not simply a question of academic interest, and to answer it is of utmost impor- 
tance. 

In Refs. [2j[3] scale-invariant solutions in d = 4 — e were presented using the two-loop 
beta functions. However, it was identified there that, in the minimal subtraction (MS) 
scheme, the solutions could disappear if the three-loop contributions to the beta functions 
were included. By performing a scheme change it was possible to make the solutions 
trustworthy at two-loop order, since the classical contribution to the beta functions does 
not participate in the scheme change. Although technically correct, this makes the examples 
found in Refs. [2,3| somewhat displeasing. 

In Ref. [3] a scale-invariant solution in d = 4 with unbounded-from-below scalar poten- 
tial was also presented. Again, the solution was obtained using the two-loop beta functions, 
and it is also subject to the criticism that the three-loop terms in the beta function may 
eliminate it. Furthermore, since the beta functions are covariant under appropriate scheme 
changes in d — 4, no scheme changes in d — 4 can modify the scale-invariant solutions, 
which can therefore disappear at three loops in any scheme. 

Thus, although scale does not necessarily imply conformal invariance in a unitary quan- 
tum field theory (QFT) with enough scalars and fermions at two loops |2], no completely 
trustworthy examples have been discovered at that order. The failure to find concrete 
examples at two loops can be understood using the results of Osborn [6j[7] and Jack & 
Osborn j8]. In Ref. [8] it is argued that, in the weak-coupling regime, RG flows are gradient 
flows at two loops. Hence, even though scale does not necessarily imply conformal invari- 
ance at two loops, the beta function monomials which could lead to concrete scale-invariant 
theories have coefficients that conspire to make all solutions conformal. Nothing forbids 
this from occurring order by order in perturbation theory. Therefore, either scale implies 
conformal invariance — and the coefficients of the beta function monomials are tightly con- 

1 For other explorations of scale without conformal invariance see Refs. 5 . 
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strained, forcing all would-be scale-invariant solutions to be conformal — or it does not — and 
recurrent behaviors exist. Either way, the answer to the original question leads to impor- 
tant implications (unexpected structure in the beta functions or the existence of recurrent 
behaviors) and the question deserves to be fully investigated. 

In this paper we compute the necessary three-loop contributions to the beta functions 
to determine if the scale-invariant solutions found in d = 4 — e are eliminated at three loops 
in the MS scheme, i.e., without an arbitrary scheme change. Our results show that the 
scale-invariant solutions are robust at three loops, and thus open the door for a d = 4 scale- 
invariant example. Indeed, since scale implies conformal invariance in pure gauge theories 
at weak coupling [TJ[8], the addition of gauge bosons in d — 4 should not qualitatively 
change the d = 4 — e results. For example, the beta function monomials exhibited below, 
which lead to an obstruction to the gradient flow interpretation of the RG flow, are not 
modified in any way by the introduction of gauge bosons. However, to fully answer the 
question in d = 4, one needs the complete three-loop beta functions of theories with matter 
and gauge fields, a computation we hope to undertake soon. 

It is important to point out that the c-theorem discussed in Refs. j^-j8], which leads to 
dc/dt = —Gufi 1 (3 J with Gjj positive-definite in the weak coupling regime, is too restrictive. 
Indeed, following Osborn [7] , the all- loop proof of the c-theorem, which implies the existence 
of a monotonically decreasing c-function which is constant only at conformal fixed points, 
must be modified once spin-one operators of dimension three are taken into account. This 
is exactly the case for non-conformal scale-invariant theories, since the virial current is such 
an operator. Taking into account the virial current, the analysis is modified as described 
in Ref. (TJ section 3], and leads to dc/dt = —(Gjj + ■ ■ ■)f3 I B J where B 1 = (3 1 — Q 1 and 
/3 1 = Q 1 for non-conformal scale-invariant theories. Thus, in its most general form the 
work of Osborn [6j[7] and Jack & Osborn [8] implies the existence of a c-function which 
is constant at conformal fixed points ((3 1 = 0) as well as on scale-invariant trajectories 
(B 1 = 0). This is in accord with our three-loop results. 

The paper is organized as follows: In section [2j we discuss the e expansion in more 
detail, showing why the scale-invariant solutions found in Refs. [20 can be destabilized at 
three loops. We then generate the most general three-loop beta function for the Yukawa 
coupling and determine which diagrams contribute to the virial current. We finally compute 
the beta function coefficients of the relevant diagrams and verify that the virial current does 
not vanish at three loops, thus demonstrating the existence of scale-invariant theories in 
d = 4 — e in a well-defined renormalization scheme. Along the way we elucidate the subtle 
effects of scheme changes in d = 4 — e dimensions. In section |3j we examine scheme changes 
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in theories with many couplings and also on scale-invariant solutions, showing that, as 
expected, physical parameters in d = 4 do not depend on the renormalization scheme. 
In section [| we elucidate the stability properties of scale-invariant solutions and explicitly 
verify that the example of Ref. [3] exhibits both attractive and repulsive directions. In 
section [5] we return to the arguments of Osborn [6j[7] and Jack & Osborn [8] and show that 
they are not in contradiction with our results. Finally, in section [6] we contrast our cyclic 
trajectories with the trajectories of Ref. (9] which were recently discussed in connection 
with the c-theorem in Ref. [l0| (see also Ref. [IT]). 



2. Establishing scale invariance 

2.1. The two-loop computation 



The results of Refs. [2j|3] were presented in an expansion in e, similar in spirit to the 
expansion that reveals the Wilson-Fisher fixed point. Let us recall here how that works. 
We consider a model with real scalar fields (f) a and Weyl spinors ipi with quartic scalar 
self-couplings \ a bcd and Yukawa couplings y a \ij- The equations for scale invariance are 

(3abcd(^,y) = Qabcd = — Qa'a^a'bcd — Qb'b^ab'cd — Qc'c^abc'd ~ Qd'd^abcd' , (2.1a) 
Pa\ij(Ky) — Pa\ij = — Qa'aVa'lij ~ Pi'iUa\i'j ~ Pj'jUa\ij' > (2-lb) 

where f3 a bcd = —d\ a bcd/dt and (3 a \ij = —dy a \ij/dt are the beta functions for the coupling 
constants^] Q a b is antisymmetric and P^ anti-Hermitian. To proceed, we solve Eqs. (2.1) 
for the coefficients of A, y, Q and P in an e expansion, 

Kbcd = ^ ^oLd 6 " Va\ij = ^ y^ij 6 * 1 '*' = ^i^ 6 "' Pi i = P if en ' 

n>l n>l n>2 n>2 

For the remainder of this section we will work with a theory of two real scalars and 
two Weyl fermions, with canonical kinetic terms and interactions described by 

V = iAi$ + ^\ 2 <f>\ + |A 3 0?02 + \^4>\<p2 + |A 5 0i02 + ikviMrti 

+ |l/202'0lV'l + ^2/301^2^2 + ^1/402^2^2 + 2/501^1^2 + 2/602^1^2 + h.C.) . (2.2) 

This is the simplest weakly-coupled unitary example in d = 4 — e with a well-behaved 
bounded-from-below scalar potential. For this model = and Q a b is 2 x 2 antisymmetric 
with Q 12 = q. 



2 With our conventions RG time increases as we flow to the IR, t — hx{pbo/y). 
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Eq. (2.1b) is to be solved first, at order e 3//2 . The result is used in Eq. (2.1a) which is 
then solved at order e 2 . This is a system of coupled nonlinear equations and, as such, it 
has many solutions y^}. and X a ]j cd , some of them consistent with unitarity and boundedness 
of the scalar potential, while others not. Additionally, some of these solutions lead to 
conformal fixed points, while others allow for nonzero q, at least in principle. 

At two-loop order solutions y^- and A^ cd of the previous order are used to solve 
Eq. (2.1b) at order e 5//2 , and Eq. (2.1a) at order e 3 . This is now a system of coupled 



linear equations] 3 from which the unknowns 



and X^abcd the unknown q 



(2) 



Va\ij alld A(2! 



,(1) 



abed are determined. For most 



is equal to zero, but for certain y^}. and A^ J cd , i.e., for possi- 



ble scale-invariant solutions, it is found be equal to a linear combination of coefficients of 
monomials in a uj. More specifically, the diagrams that contribute to q at two loops are 
shown in Fig. [TJ 




Fig. 1: Diagrams that contribute to q at two-loop order. 



Let b± and bi be their coefficients in the beta function for the Yukawa coupling, 

(167r 2 ) 2 /3^ loop) D hy^y^y^jXabcd + b 2 y b \ i:j X bcde X acde . 

Then (we omit the prefactor here since it is not relevant for the discussion), 

g (2) oc 61 + 246 2 , (2.3) 

and because bi = —2 and 62 = we find q^ = 0. As we already mentioned, the failure to 
find trustworthy non-conformal scale-invariant solutions at two loops can be explained by 
the gradient flow property of the RG flow at weak coupling described in Ref. [8). Note that 
here, contrary to the case of conformal fixed points, 7^ at two-loop order. However, 
the three- loop contributions to the beta functions can very well conspire to set g'- 3 -' = 0, 

3 For all higher orders in e one only gets systems of coupled linear equations. 
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and thus restore conformal invariance. (As we will demonstrate in the next subsection, this 
does not happen. The fact that = is merely an accident.) 

An interesting observation at this point is that if = were not an accident, then 
that would directly imply that the conformal symmetry somehow relates vertex corrections 
and wavefunction renormalizations. This is obvious from the fact that the first diagram 
in Fig. [I] contributes to the residue of the 1/e pole of Z y , while the second to the residue 
of the 1/e pole of Z^. This would be reminiscent, e.g., of the Ward identity for charge 
conservation in QED. 

Now, in d = 4 — e we can perform a scheme change and get q^ ^ 0, for, in d = 4 — e, 



the equations for scale invariance (2.1) are not covariant under scheme changes. The reason 
is the classical origin of the first term in the beta functions, which does not take part in 
the scheme change. Consequently, the d = 4 — e results of Refs. [2j[3] (after the scheme 
change) suggest that the scale-invariant solutions can be trusted, solely, however, due to 



the unavoidable noncovariance feature of Eqs. (2.1) in 4 — e dimensions. On the contrary, 
in d — 4 there is no classical contribution to the beta functions, and so the results of 
Ref. (3) in d = 4 are not trustworthy in any scheme, i.e., they are not robust against the 



inclusion of the three-loop terms in the beta functions. (For further discussion of the effects 



of scheme changes see subsection 3.3 ) 



2.2. The three-loop computation 

There are a few motivations to go to the next order in e. Firstly, the scheme change that 
establishes scale-invariant solutions in d = 4 — e is arguably undesirable, and makes the 
phenomenon look like an artifact of 4 — e. It is critical to see the effect in a well-defined 
scheme, which for example should have a well-behaved e — > 1 limit and thus might be of 
relevance in strongly-coupled d = 3 models. Secondly, and more importantly, including 
the next order makes the extension to d — 4 possible. Indeed, in d — 4 the role of 
e is played by the critical value of the gauge coupling at the Banks-Zaks fixed point, 
and the corresponding terms in the beta functions do contribute to the scheme change. 
Consequently, the scheme change that reveals scale without conformal invariance in d = 4— e 



cannot help in d = 4, and the next order is needed. As in Ref. 12 we will perform our 
computation in dimensional regularization with the MS scheme. 

There is a large number of diagrams that contribute to (3 a \ij at three loops. (We use 
the Mathematica package FeynArts to automatically generate all required diagrams.) Each 
diagram corresponds to a unique combination of coupling constants, and an examination 
of all of them reveals that those that contribute to q^ are the ones shown in Fig. [2J 
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Fig. 2: Diagrams that contribute to q at three-loop order. 
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Note that the diagrams in Fig. [2] are specific to the example of this subsection. More 
complicated models might involve more diagrams. It is interesting to point out that very 
few of the ~ 200 diagrams in f3 a uj contribute to g^, and that the ones that do involve 
both Yukawa and quartic vertices. The same situation is encountered at two loops, and 
we conjecture that it holds to all orders in perturbation theory. This is also motivated by 



comments in Ref. 13 regarding the "interference" between successive loop orders in the 



calculation of a potential for a gradient flow (see also section [5] below). 

The diagrams of Fig. [2] have simple poles in e and so they contribute to the Yukawa 
beta function at three loops: 

mnVe\nj^abde + ' ' ' + Cl2yb\ij^bcde^cdfg^aefg- 



(2.4) 



The three-loop analog of Eq. (2.3) is then (again omitting the prefactor) 

g (3) oc -71 + 3(ci + 2c 2 + 2c 3 + c 4 + 2c 5 + 4c 6 + 8c 7 ) + 4(c 8 + 2c 9 + 3c w + 4cn + 58ci 2 ) 

where the constant piece comes from contributions to from the previous order. 

To compute these three-loop diagrams we implemented the algorithm of Ref. 
There, IR divergences are regulated by introducing a common mass parameter to all propa- 
gator denominators, and the calculation proceeds with properly choosing a loop momentum, 
regarding it as large, and expanding with respect to it the remaining two-loop subintegral, 
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for which the chosen momentum is external. Remarkably, the authors of Ref. 14 manage 



to construct explicit formulas for the pole parts of all three-loop scalar integrals. The im- 
plementation of their algorithm is straightforward, but one must be very careful to take 
into account all required counterterms, including the ones introduced by the IR regulator. 
To test our implementation, we verified the two- loop result of Ref. [l2j for (3 a \ij, and also 
part of the three-loop result for the beta function of the quartic coupling in a multi-flavor 
theory of scalars found in Ref. j8]. We also performed explicit computations of a couple of 
diagrams. 

From the diagrams of Fig. [2] we find 



Cl 
C7 



i 

2' 



C-2 

C 8 



3 
2 • 



C 3 

c 9 



1 

2' 



ClO 



c 5 

Cll 



1 

2' 



_5_ 

32' 



06 
Cl2 



3 
2 • 



16' 



Restoring the prefactor, then, Eq. (2.4) gives 



,(3) 



^(757-3^52 953) 



108 300 



7 x 10 -5 . 



l We would like to thank M. Misiak for pointing us to this reference. 
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Since ^ no scheme change is needed here to establish the existence of theories that 
are scale but not conformally invariant. Because of this we expect that theories in d — 4 
can also display scale without conformal invariance. 

To summarize, it is important to emphasize that the distinction between scale-invariant 



and conformal solutions of Eqs. (2.1) at the two-loop level is that, for the latter, q^- 3 ^ = 
already at two loops. Higher loops are expected to slightly modify the critical values of 
the couplings, while preserving q = 0. But there are solutions for which g(- 3 ) already 
at two loops. As a result, the nature of these solutions is uncertain, and a higher-loop 
calculation is needed. Even without that calculation, though, it should be clear that not 



all solutions to Eqs. (2.1) can be declared conformal with the same confidence. There are 
doubts for some of them, and in d = 4 — e the change of scheme done in Refs. [2j|3], although 
undesirable, suggests that they are indeed not conformal. The three-loop computation we 
present here shows that there is no need for a scheme change: non-conformal scale-invariant 
solutions exist in a "standard" scheme. This offers strong indications that they also exist 
in d — 4. 



3. Renormalization-scheme changes 

3.1. Scheme changes and conformal fixed points: the one-coupling case 

Let us first review the effects of scheme changes in conformal theories. The simple case 
of a theory with only one coupling has been investigated long ago in Ref. [15]. Under a 
scheme change, the coupling g and the wavefunction renormalization Z(g) become 

g^g(g)=g + 0(g 3 ), 
Z^(g)^Z^(~9) = Z^(g)F(g), ^ 

where F(g) = 1 + 0(g 2 ) and F^O for all g. In the new scheme g is equal to g at lowest 
order since the coupling is unambiguous at the classical level. The same is true for the 
wavefunction renormalization as well. Therefore, since^] 



7(9) = -z-/ 2(9) ^ 



5 We use B = Z 1/2 (ff)</r 



the new beta function and anomalous dimension are related to the old beta function and 
anomalous dimension through 

9 BF( \ (3 - 2) 

^g) =1 (g) + F~\g)P{g)-^. 

Although the RG functions depend strongly on the renormalization scheme, properties that 
have physical consequences must be independent of the scheme. Such properties are: 

(I) The existence of a conformal fixed point; 

(II) The anomalous dimension at a conformal fixed point, which determines the scaling 
behavior of Green functions; 

(III) The first derivative of the beta function at a conformal fixed point, which determines 
the sigri0 and rate of approach of the coupling to the conformal fixed point and thus 
modifies asymptotic formulae; 

(IV) The first two coefficients in the beta function, which govern the UV or IR asymptotics 
of the coupling; 

(V) The first coefficient in the anomalous dimension, which controls the scale factor of 
the field in the far UV or IR. 



These properties all follow from Eqs. (3.2) and the form of g(g) and F(g). 



3.2. Scheme changes and conformal fixed points: the multi- coupling case 

When the theory has more than one coupling, a scheme change still transforms the coupling 
vectoiQ g 1 and the wavefunction renormalization matrix Z/(g) as in (3.1) but, due to the 



vector and matrix character of the coupling and wavefunction renormalization respectively, 
the new wavefunction renormalization is modified by a matrix F T J (g) through 

Z^{g) -+ Z^(g) = Z 1/2 (g)F(g). 



6 Note that the sign determines the character (attractive or repulsive) of the conformal fixed point. 
7 Capitalized indices run through all couplings. For matrices we use, e.g., Q/ for both Q a b and Py. 
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Thus, under a scheme change, one has 

dg 1 



T (9) = P J (g) 



dg v 

[F-\gUg)F{g)\; + 



F-\g)P K {g) 



dF{g) 



(3.3a) 
(3.3b) 



It is easy to see that, in the multi-coupling case, properties (I) and |(V) are still scheme 



independent. Property (II) is of course modified so that only tr7 and det7, and so the 



eigenvalues of 7, are scheme- independent. Property (III) is also modified since 

dp J (~g) 



dg 



hi 



dg K dp L {g) d~g J + dg K ^ d 
dg 1 dg K dg L dg 1 dg L 



dg 
dg 



7.1 



A 



(3.4) 



such that at a conformal fixed point the eigenvalues of dP J (g)/dg I are independent of the 
scheme. This is expected because d(3 J /dg 1 = 7 7 J , where 7/ is the anomalous-dimension 



matrix of the operators sourced by the appropriate couplings. Therefore, Eq. (3.4) can be 
seen as an extension of Eq. (3.3b) with F = dg/dg. 



Finally, if the one-loop beta function for one coupling depends on other couplings, 



property (IV) is no longer true [3j| — only the first coefficient in the beta function is scheme- 



independent, although the UV or IR asymptotics of the couplings are the same in any 
scheme. 

3.3. Natural scheme changes and scale-invariant trajectories 

It is interesting to see how scale-invariant solutions behave under scheme changes]^] Here we 
will distinguish between two types of scheme changes, which we dub natural and unnatural. 
A natural scheme change transforms the couplings as 

^abcd ^abcd = ^abcd + Vabcd ) 

Ua\ij ~^ Va\ij = Va\ij + ^a\ijj (3-5) 
Va\ij ^ Va\ij ^a|«j ^a\iji 

such that all couplings transform covariantly with respect to the symmetry group of the ki- 
netic terms. MS and variants are examples of this — it occurs, e.g., every time one dresses a 
Feynman diagram topology with couplings. Unnatural scheme changes spoil the covariance 
of equations. 



The discussion of this subsection applies to scheme changes under which Eqs. (2.1 ) transform covariantly. 



Since the analysis for gauge fields is straightforward, gauge fields are omitted for simplicity. 
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We can now show that entries of Q and P, which determine, e.g., the frequency on a 
cyclic trajectory, are scheme-independent for natural scheme changes. Indeed, if the scheme 
change is natural, then the time evolution of i] and £ on a scale-invariant trajectory is given 
by 

Vabcdit) = (e Qt ) a , a (e Qt )v b (e Qt ) c , c (e Qt ) d/d ri alb , c , d ,(0), 

(3.6) 

6.|tf(t) = (e Q *)a'a(e P *)«(e pt ),v i e a , K , i ,(0), 



and so 



= Qa'aVa'bcd + permutations, = Q a > a £,a'\ij + Pi'i€a\i'j + Pj'jU\ij'- (3-7) 



On a scale-invariant trajectory Eqs. (3.5) give 



5 _ n drjabcd n _ ^ dja\ij / q m 

Pabcd — ^dabcd ^ — , Pa\ij ~ ' a\ij 



and we can use Eqs. (3.7) to obtain 



Pabcd = -Qa'aK'bcd + permutations, 

Pa\ij Q a' aVa'\ij Pi'il/a\i'j f j' 1 'j% la\ij' '• 

Hence, Q and P are scheme-independent for natural scheme changes. 

As a result of our analysis the existence of scale-invariant trajectories does not depend 



on the renormalization scheme. As expected, then, property (I) is easily extended to include 
non-conformal scale-invariant trajectories. 

Focusing on scalar anomalous dimensions (the argument can be easily repeated for 
fermion anomalous dimensions), property (II)| can also be generalized to scale-invariant 



theories. Indeed, for natural scheme changes on a scale-invariant trajectory Eq. (3.3b) 
becomes 

laM = [F- 1 (9h(g)F(g)] ab +{F-\g)[Q,F(g)}} ab 
since —dF(g)/dt = [Q,F(g)]. One can then immediately see that (using matrix notation) 

m + Q = F'\g)h(g) + QW(g), 

so that the eigenvalues of 7 + Q are scheme-independent. This is in accord with expecta- 
tions: in Ref. [4] it was shown that the behavior of two-point functions is determined by 
the eigenvalues of 7 + Q, which are therefore expected to be scheme-independent. 
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Since property (II) can be generalized to scale-invariant theories, the same is expected 
for property (III) due to d/3 J /dg 1 = j/. Indeed, Eq. (3.4) becomes 

d(3 J 



)7,I 



F-\g)f g F{g) 



+ {F- 1 (g)[Q,F(g)}} I J , 



where F = dg/dg, which gives (again using matrix notation) 

'dp 



dg 



+ Q 



F(g)- 



Therefore, the eigenvalues of d(3/dg + Q = d({3 — Q)/dg (since Q = —gQ) are scheme- 
independent. It is interesting to note that the eigenvalues of 8(3/ dg + Q are expected 
to determine the character (attractive, repulsive, etc.) of scale-invariant trajectories, and 
so one of them should be zero — that is indeed the eigenvalue corresponding to the (left) 
eigenvector (3 1 = Q 1 . This is because /3 7 = Q 1 generates a motion along the scale-invariant 
trajectory, not away from it, as can be seen directly from 



1 



d/3 



dg 



7 + Qi J 



dQ 1 
dt 



+ Q'Q/ = 0. 



Finally, properties (IV) and |(V) in the multi-coupling case are trivially extended to 
scale- invariant theories since they do not rely on the existence of scale- invariant trajectories 
(or conformal fixed points). 

To summarize, the scheme-independent properties (I-V) can be generalized to: 

(I') The existence of conformal fixed points and scale-invariant trajectories; 

(II') The eigenvalues of 7 + Q at conformal fixed points and scale-invariant trajectories; 

(III 7 ) The eigenvalues of d/3/dg+Q at conformal fixed points and scale-invariant trajectories; 

(IV) The first coefficient in the beta functions; 

(V) The first coefficient in the anomalous-dimension matrix. 



4. Stability properties 

4-1- General discussion 

It is of interest to study the stability of scale-invariant solutions under small deformations. 
Such an analysis determines the character of a particular scale-invariant solution, which 
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can have (IR) attractive and/or repulsive deformations. In this section we will describe the 
properties of all possible scale-invariant solutions. The corresponding results for conformal 
fixed points are recovered by setting Q = in the equations below. To simplify the 
equations, matrix notation is used throughout this section. 

Since non-conformal scale-invariant solutions exhibit non-trivial RG flows, it is natural 
to disentangle the two contributions to the flow of the deformations, i.e., the expected 
contribution from the non-conformal scale-invariant solution, and the actual contribution 
from the deformations which we want to analyze. The appropriate quantity to study 
is thus Sg(t) = [g(t) — g*(t)]e~ Qt , where = g*(0)e Qt is a scale-invariant solution, 

fi\g= g *(t) = Q(t)- The quantity 5g{t) determines the behavior of the deformations as a 
function of RG time in a "comoving frame", i.e., modulo the expected non-conformal scale- 
invariant solution RG flow. Note that, although for non-conformal scale- invariant solutions 
the choice of g*(0) in Sg(t) = g(t)e~^' — g*(0) is arbitrary]^] in order to study the behavior 
of small deformations one should first fix a g*(0). 

To proceed further it is necessary to Taylor expand the beta functions around the 
appropriate scale- invariant solution 



m = P\ a=9f{t) + W)-g,{t)]j- g 



9=9*(t) 



+ ... = Q(t)+Sg(t)^- 



e Qt + ... 



9=9*(0) 



where the last equality follows since —d(d(3 /dg)/dt = [Q,d(3/dg] on the scale-invariant 
solution. Note that in order to disentangle the two contributions to the flow, the above 
Taylor expansion is RG-time dependent. It is now straightforward to write down, at lowest 
non-trivial order, the system of (linear) differential equations that the deformations must 
satisfy: 

= W) - Q(t)}e~ Qt + 5g(t)Q = 6g(t)S + ■■■, (4.1) 



dt 
where 



+ Q 

g=g*(o) 



S = 

is the stability matrix. It is obvious that 5g (t) is the appropriate choice of variable that 



allows a separation of the RG flow contributions, for all RG-time dependence in Eq. (4.1) 



comes solely from 5g(t). Note, moreover, that Eq. (4.1) implies that the behavior of the 
deformations 5g{t) is dictated by the eigenvalues of S which, as we showed in the previous 



9 Any two points on a non-conformal scale-invariant trajectory are physically equivalent due to scale 
invariance. 



13 



section, are scheme-independent (property (III'))- The solution to the system of differential 
equations (4.1) is simply 

6g(t) = 6g(0)e- st + • • • (4.2) 

and one can easily see that positive (respectively, negative) eigenvalues of the stability 
matrix S correspond to IR attractive (respectively, repulsive) deformations. As usual, the 
fate of deformations related to vanishing eigenvalues cannot be determined from Eq. (4.2)— 



for vanishing eigenvalues it is necessary to go to higher order in the Taylor expansion (4.1). 
However, as already mentioned, non-conformal scale-invariant solutions exhibit one special 
(left) eigenvector Sg(0) oc Q(0) with vanishing eigenvalue which represents a deformation 
along the scale-invariant solution. For this special deformation the full solution Sg(t) = 
[g*(t ± St) - g*(t)]e~ Qt = g*(Q)[e ±QSt - 1] = TQ(0)5t + ••• is RG-time independent as 
expected, since it corresponds to a flow along the RG scale-invariant trajectory. 

The previous analysis is a generalization of the similar analysis done for conformal 
solutions where Q — 0. Note that the special (left) eigenvector Sg(0) oc <2(0) does not 
exist for conformal fixed points, as expected since conformal solutions do not exhibit any 
non-trivial RG flow. 



4-2. The example 

We can now use the results discussed above to investigate the behavior of small deforma- 
tions away from scale-invariant solutions. To this end it is natural to use an e expansion 
for the stability matrix S and its eigenvalues x m , 

S = J2S^el x m = Y J ^ ) ^- 

n>2 n>2 

The form of the expansion is dictated by the form of the beta functions in the stability 
matrix. 

The eigenvalues of the stability matrix are the roots of the characteristic polynomial 
det(xl — S) which can also be expanded in e. To lowest order the characteristic polynomial 
simplifies and the eigenvalues are solutions of 

det(x (1) l - S w ) = 0. (4.3) 

Since there are only seven non-vanishing independent couplings (Xi > ,..,5,yi ; 2 i n ( ]2-2[ )) at 
a generic point on the non-conformal scale-invariant solution described in section [2] (see 
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Ref. [3] for details), Eq. (4.3) for the corresponding couplings is 



7(7 — -\ \ [ 7 5_ V52953 4 , 1894+V52953 3 _ 240768-335V52953 2 
Z ^ ^ ^ 57 Z " r 475 ^ 135375 Z 

421203- 1573V/52953 , 136(757y / 52953- 158859) \ = n 
225625 64303125 J 

which cannot be solved by factorization into radicals. (To avoid clutter we define z = x^\) 
A numerical solution gives five positive, one negative and one vanishing eigenvalue: 

Z&2A, z = l, z«0.99, 2^0.74, z ns 0.095, z « -0.19, z = 0. 

The positive eigenvalues show that the scale-invariant solution is IR attractive in several 
directions. We thus expect that the limit cycle can be reached by an appropriate defor- 
mation of a theory defined at a UV conformal fixed point, although, to be certain, a more 
thorough analysis is necessary. 



5. On the proof of the c-theorem at weak coupling 

As discussed in the introduction, our three-loop results do not contradict the work of 
Osborn (6|[7j and Jack & Osborn |8|. Focusing on Ref. [7], Osborn proved that RG flows 
are gradient flows at two loops in the weak coupling regime. Lifting the theory to curved 
space with spacetime-dependent couplings, Osborn showed that Weyl consistency conditions 
lead to 

with Gjj positive-definite in the weak coupling regime, thus forbidding the existence of 
recurrent behaviors at all loops. From the general analysis of Ref. [7] it would thus seem 
that scale-invariant trajectories are forbidden to all orders in perturbation theory. However, 
the analysis of Ref. [7] leading to Eq. (5.1) is too restrictive — it does not allow for spin- 



one operators of dimension three, i.e., it does not include the possibility of non-conformal 
scale-invariant theories. 

The more general analysis, also performed by Osborn in Ref. [7], includes possible 
spin-one operators of dimension three, which are related to the symmetry group of the 
kinetic terms. Such an analysis is done by promoting the related symmetry of the kinetic 
terms — for example the symmetry of the kinetic terms generated by the virial current, the 
natural spin-one operator of dimension three for scale-invariant theories — to a symmetry 
of the interacting theory. This is implemented by allowing the couplings to transform 
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appropriately under a change generated by the spin-one operators of dimension three and 
by introducing background gauge fields to render the symmetry local. Then, assuming that 
the regularization procedure preserves local gauge invariance, Osborn's Weyl consistency 
conditions and current conservation show that 

| = _^|l = _ (G „ + .. W (5.2) 

where B 1 = (3 1 — Q 1 . Note that B 1 = is precisely the condition for scale invariance. Thus, 
by allowing non-conformal scale-invariant theories from the start, the work of Refs. pj-p] 
implies the existence of a c-function whose RG-time derivative vanishes at conformal fixed 
points as well as on scale-invariant trajectories. Note, moreover, that the c-function might 
not be monotonically decreasing due to the extra contributions to dc/dt represented by the 
ellipsis in Eq. (5.2). 

Note that, by promoting the symmetry of the spin-one operators of dimension three to 
a symmetry of the interacting theory, it is natural to demand regularization and renormal- 
ization schemes that satisfy the newly promoted symmetry. This also explains the special 
status of the natural renormalization schemes defined in the previous section. 

Finally, it is interesting to see why the interference between quartic coupling one-loop 
beta functions and Yukawa coupling two-loop beta functions proposed by Wallace & Zia [13) 
as a possible obstruction to the gradient flow interpretation of the RG flow is circumvented 
by the introduction of the metric. Focusing on the problematic monomials in a possible 
c-function, 

c D di tr(y* a y b y* c y d )X abcd + d 2 tr(y* a y b )X acde X bcde , 
the related contributions to the beta functions at one and two loops respectively are 
dc 



dXabcd 

dc 
dy a 



D di tr(y* a y b y*y d ) + 2d 2 tr(y d y e )\ abce + permutations, 
D 2d 1 y b y*y d X abcd + d 2 ybX ac deXbcde- 



Comparing with the true beta functions, 

1 11 

(3ab!d° P) 3 -^^tehfaVbVlVd) + -tr(^y J A a6ce + permutations, 

(2 1) ^ 11 

P a P D — Jy^^ybVcVd^abcd + (167T 2 ) 2 12 ^ acde ^ 6c£fe ' 

it is straightforward to see that the metric can account for the loop mismatch since d 2 jd\ = 
— 1/12 for both beta functions, as pointed out in Ref. (8l. Note that the conditions for a 
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gradient flow interpretation of the RG flow introduced at higher orders are ever more 



constraining due to the large number of diagrams and it is plausible that they are not 



satisfied, as our three-loop computation shows. The interference argument of Wallace & 



Zia 13 prevails at three loops, although for a complete investigation the knowledge of the 
full three-loop beta functions is necessary. Interestingly, the interference between the (n— 1)- 
loop quartic-coupling beta function and the n-loop Yukawa beta function also explains why 
the n-loop quartic-coupling beta function is not necessary to argue for the existence of 
scale-invariant theories at n-th order in perturbation theory. 



6. Cyclic trajectories and the c-theorem 

It is important to note that the existence of recurrent behaviors in RG flows in d — 4 
does not contradict all versions of the c-theoremj^] In particular, the weak version of the 
c-theorem, where two conformal fixed points connected by an RG flow satisfy the inequality 

ouv - «m > (6.1) 

with a the conformal anomaly (see, for example, Ref. |17|){^]is consistent with scale without 
conformal invariance. Even the stronger version of the c-theorem, where there exists a 
local function which is monotonically decreasing along non-trivial RG flows, is compatible 
with recurrent behaviors as long as the c-function is constant on scale-invariant trajectories. 
Only the strongest version of the c-theorem is violated by the existence of limit cycles and 
ergodicity; a gradient flow interpretation of RG flows is impossible for theories in which 
scale does not imply conformal invariance. 

Since theories exhibiting limit cycles or ergodicity are scale-invariant, it is reasonable 
to expect the interpolating c-function to be constant on scale-invariant trajectories. Any 
such interpolating function is invariant under the symmetry group of the kinetic terms, 
i.e., it does not carry scalar or fermion indices. Thus, in a natural scheme, all the explicit 
RG-time dependence disappears on a scale-invariant trajectory. This is the behavior that 
is intuitively expected of the c-function, which should be some measure of the number of 
massless degrees of freedom of the theory. Therefore it must be constant on scale-invariant 
trajectories since any two points on such trajectories are physically equivalent. 



This was already noticed in Ref. 



1G 



10 

11 For a more extensive discussion see Ref. jij. 

12 A claim for the proof of the inequality (|6.1|) appeared recently in Ref. 
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(see also Ref. 



19 
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This behavior is very different from that encountered on cyclic flows described in Ref. [9] 
and recently discussed in association with the c-theorem in Ref. [To] (see also Ref. [IT]). 
In Ref. 1 10 1, the authors argue that monotonic RG flows can be simultaneously cyclic if 
one allows for a multi-valued interpolating c-function. This is fundamentally different from 
recurrent behavior with continuous scale invariance. As mentioned above, the interpolating 
c-function must be constant on scale-invariant trajectories. Moreover, the examples cited in 



Ref. 10 exhibit one feature, turning points, which does not appear on continuously scale- 
invariant trajectories. Turning points are peculiar locations in coupling space: the beta 
functions vanish there, but the first derivative of the beta functions diverges. Consequently, 
RG flows can overshoot turning points. In contrast, all existing continuously scale-invariant 
examples are well-defined smooth weakly-coupled theories, and thus do not display turning 
points. The existence of turning points on cyclic flows is a reflection of the possibility of 
multi-valued c-functions which are monotonically decreasing along the flow. Here we want 



to stress that the physics of cyclic flows with turning points as described in Ref. 10 is very 



different from that of recurrent behaviors with continuous scale invariance. It is therefore 
very unlikely that monotonically decreasing multi- valued c-functions exist on scale-invariant 



recurrent behaviors as suggested in Ref. 10 



7. Conclusion 

Does scale imply conformal invariance in unitary relativistic QFTs? The answer is negative 
in d = 4 — e. Although a similarly conclusive statement in the d = 4 case cannot yet be 
made, we strongly believe that the answer there is also negative. There are no physical 
arguments on which one can rely to forbid non-conformal scale-invariant theories. Instead, 
one simply needs to compute the beta functions and explore the different regions in coupling 
space. That an example of a scale- invariant theory which is not conformal eluded the 
physics community for so long is easily explained by the complexity of the problem: to see 
non-conformal scale-invariant theories, one must go to three loops, and the beta functions 
at three loops in the most general QFT are not known. 
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